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Abstract. We consider a non relativistic charged particle in a ID infinite square potential well. 
This quantum system is subjected to a control, which is a uniform (in space) time depending electric 
field. It is represented by a complex probability amplitude solution of a Schrodinger equation on a 
ID bounded domain, with Dirichlet boundary conditions. We prove the almost global approximate 
stabilization of the eigenstates by explicit feedback laws. 
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1. Introduction. 

1.1. Main result. As in [30l [6l [8] . we consider a non-relativist charged particle 
in a one dimensional space, with a potential V{x), in a uniform electric field t 
u(t) € M. Under the dipole moment approximation assumption, and after appropriate 
changes of scales, the evolution of the particle's wave function is given by the following 
Schrodinger equation 

dm a^* 

i — {t, x) ^- — {t,x) + {V{x) - u{t)x)m{t, x). 

We study the case of an infinite square potential well: V{x) = for x G / := 
(—1/2, 1/2) and V{x) = H-oo for x outside /. Therefore our system is 



i — {t,x) = -- — {t,x)-u{t)x-^{t,x), xGl (1.1) 



m{0,x)^mo{x), (1.2) 
«'(t,±l/2) = 0. (1.3) 



It is a nonlinear control system, denoted by (S), in which 

• the control is the electric field u : M+ M, 

• the state is the wave function : K_|_ x / ^ C with \E'(i) e § for every t ^ 0, 
where § := {(ys e L'^{I;C)]\\ip\\L2 ^ 1}. 

Let us introduce the operator A defined by 

D{A):=iH'nH^){I,C), A^:^-\^, 
and for s € M the spaces 

i/fo)(/,C) :^Z?(A^/2). 
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The following Proposition recalls classical existence and uniqueness results for the 
solutions of p.ip - (|1.2p - (|1.3p . For sake of completeness, a proof of this proposition is 
given in the Appendix. 

Proposition 1.1. Let *o e §, T > anrf u e C°([0,r],R). There ex- 
ists a unique weak solution of il.l\) - n^) - U.3\) . i.e. a function '5 e C°([0,r],§) fl 
Ci([0,T],ff(-^(/,C)) .such that 

-^{t) ^ e-'^^'^o + i [ e-^'^'^^-'^u{.s)x'<i'{s)ds in L^{I,C) for every te[0,T]. (1.4) 



and then m.l\) holds in H^^^ (/, C) for every t G [0, T] . 

If, moreover, VI/q S {H^ H Hq)(I,C), then is a strong solution i.e. ^' G 
C°([0,r],(i/2nFi)(J,C))nCi([0,r],L2(7,C)), the equality {flp holds in 1^(1, C), 
for every t € [0,T], the equality i 1.2(1 holds in C\ Hq{I,C) and the equality il.3]} 
holds for every t € [0, T] . 

The weak (resp. strong) solutions are continuous with respect to initial conditions 
for the C°([0, T],L^) -topology (resp. for the C°([0, T],H^n H^) -topology.) 

The symbol (., .) denotes the usual Hermitian product of L'^{I,C) i.e. 



{V-.Cl j ^Lp{x)S,{x)dx. 

For (7 e M, we introduce the operator A„ defined by 

D{A,) -.^ n Hl){TX), A,^ := - ax^. 

It is well known that there exists an orthonormal basis ((/)fc.cr)fceN' of L^(/,C) of 
eigenvectors of A„ : 



eH^n H^{I, C), A^(j)k^^ = \k,a(t>k,a 



where {\k,c!)k&i* is a- non-decreasing sequence of real numbers. For s > and a S 
we define 

equipped with the norm 

/ oo \ 1/2 



\k=l J 



For kGW and ct e M, we define 

Ck,a ■.= {q^k,ae'';e e [0,27r)}. 

In order to simplify the notations, we will write (jfk, Afc, Ck instead of (l)k,Q, A^.o, Ckfl- 
We have 

^ _ fc^TT^ 0. — / V2cos(fc7rx), when k is odd, ,^ 
2 ' I V^sm{knx), when k is even. 

The goal of this paper is the study of the stabilization of the system (S) around the 
eigenstates (pk.a- More precisely, for fc e N* and cr e M small, we state feedback laws 
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u = Ufe,CT(*) for which the solution of prTjl - pr^ -pTH )) with u{t) = Uk,a{'^{t)) is such 
that 

Hm sup distL2 (j^q {^{t), Ck,a) 

t — >+oo 

is arbitrarily small. We consider the convergence toward the circle Ck,a because the 
wave function 5" is defined up to a phase factor. For simplicity sakes, we will only 
work with the ground state 4>i^a- However, the whole arguments remain valid for the 
general case. 

Note that, even though the feedback stabilization of a quantum system necessi- 
tates more complicated models taking into account the measurement backaction on 
the system (see e.g. [201 dH [2S] ) , the kind of strategy considered in this paper can be 
helpful for the open-loop control of closed quantum systems. Indeed, one can apply 
the stabilization techniques for the Schrodinger equation in simulation and retrieve 
the control signal that will be then applied in open-loop on the real physical system. 
As it will be detailed below, in the bibliographic overview, such kind of strategy has 
been widely used in the context of finite dimensional quantum systems. 

The main result of this article is the following one. 

Theorem 1.2. Lei r > 0, s > 0, e > 0, 7 G (0, 1). There exists a** ^ a** {T, s) > 
such that, for every a e (— cr**, cr**), there exists a feedback law Va-,r,s,e,'yi^) such 
that, for every *o e S n (iJ^ n n Hf^^){I, C) with 

ll^olkf ^ r and |(*o, 01,^)1 > 7, 

the Cauchy problem fTTp-fLS^-fL^^ 'with u{t) — (J WCT,r,s,e,7(^') has a unique strong 
solution, moreover, this solution satisfies 

limsup dist]^-2[^i(t),Ci_a) ^ £■ 

t — ^+C30 

For cr 7^ 0, the feedback law will be given explicitly. For a — the feedback law 
will be given by an implicit formula. The Theorem 11.21 provides almost global approx- 
imate stabilization. Indeed, any initial condition ^Pq G § such that 'I'o € C) for 
some s > and (4'o, 4>i,a) 7^ can be moved approximately to the circle Ci^o-, thanks 
to an appropriate feedback law. We will see that the assumption "^fQ G H'' {!,€.), for 
some s > 0"is not necessary for doing that. In fact, even for a 'So only belonging to 
§, we can find the appropriate feedback law as a function of the initial state ^Pg. 

Notice that, physically, the assumption {'i'o, (pi ^a-) 7^ is not really restrictive. 
Indeed, if {^0, = 0, a control field in resonance with the natural frequencies of the 
system (the difference between the eigenvalues corresponding to an eigenstate whose 
population in the initial state is non-zero and the ground state) will, instantaneously, 
ensure a non-zero population of the ground state in the wavefunction. Then, one can 
just apply the feedback law of the Theorem 11.21 

1.2. A brief bibliography. The controllability of a finite dimensional quantum 
system, tjj^" = {Hq + u{t) iJi)* where e and Ha and Hi are N x N Her- 
mitian matrices with coefficients in C has been very well explored [33l ESI HI El [35] . 
However, this does not guarantee the simplicity of the trajectory generation. Very 
often the chemists formulate the task of the open-loop control as a cost functional 
to be minimized. Optimal control techniques (see e.g., [31j ) and iterative stochastic 
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techniques (e.g, genetic algorithms [23| ) are then two classes of approaches which are 
most commonly used for this task. 

When some non-degeneracy assumptions concerning the linearized system are 
satisfied, j26j provides another method based on Lyapunov techniques for generating 
trajectories. The relevance of such a method for the control of chemical models has 
been studied in ^7\. As mentioned above, the closed- loop system is simulated and 
the retrieved control signal is applied in open-loop. Such kind of strategy has already 
been applied widely in this framework [9l [32] ■ 

The situation is much more difficult when we consider an infinite dimensional 
configuration and less results are available. However, the controllability of the sys- 
tem (ll.ll) - (|1.2p - (|1.3p is now well understood. In 34], the author states some non- 
controllability results for general Schrodinger systems. These results apply in partic- 
ular to the system p.ip - (|1.2p - (|1.3p . However, this negative result is due to the choice 
of the functional space that does not allow the controllability. Indeed, if we consider 
different functional spaces, one can get positive controllability results. In [6], the 
local controllability of the system ()l.ip - ()1.2|) - (|1.3p around the ground state ^i,(t, for 
(T small is proved. The case cr ^ is easier because the linearized system around (pi^a- 
for a small is controllable; this case is treated with the moment theory and a 
Nash-Moser implicit functions theorem. As it has been discussed in [30j, the case of 
cr = is degenerate: the linearized system around (j>i is not controllable. Therefore, 
in this case, one needs to apply other tools, namely the return method (introduced 
in [11]) and the quantum adiabatic theory [5. In [8], the steady-state controllability 
of this nonlinear system is proved (i.e. the particle can be moved in finite time from 
an eigenstate (f)k to another one (j^j)- The proof relies on many local controllability 
results (proved with the previous strategies) together with a compactness argument. 

Concerning the trajectory generation problem for infinite dimensional systems 
still much less results are available. The very few existing literature is mostly based 
on the use of the optimal control techniques [H [5] . The simplicity of the feedback law 
found by the Lyapunov techniques in [26[ [7] suggests the use of the same approach for 
infinite dimensional configurations. However, an extension of the convergence analysis 
to the PDE configuration is not at all a trivial problem. Indeed, it requires the pre- 
compactness of the closed-loop trajectories, a property that is difficult to prove in 
infinite dimension. This strategy is used, for exemple in [14j . 

In [24j, one of the authors proposes a Lyapunov-based method to approximately 
stabilize a particle in a 3D finite potential well under some restrictive assumptions. 
The author assumes that the system is initialized in the finite dimensional discrete 
part of the spectrum. Then, the idea consists in proposing a Lyapunov function 
which encodes both the distance with respect to the target state and the necessity of 
remaining in the discrete part of the spectrum. In this way, he prevents the possibility 
of the "mass lost phenomenon" at infinity. Finally, applying some dispersive estimates 
of Strichartz type, he ensures the approximate stabilization of an arbitrary eigenstate 
in the discrete part of the spectrum. 

Finally, let us mention other strategies for proving the stabilization of control 
systems. One can try to build a feedback law for which one has a strict Lyapunov 
function. This strategy is used, for example, for hyperbolic systems of conservation 
laws in [15] , for the 2-D incompressible Euler equation in a simply connected domain 
in [12], see also [17] for the multiconnected case. For systems having a non control- 
lable linearized system around the equilibrium considered, the return method often 
provides good results, see for example [llj for controllable systems without drift and 
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[TS] ) for Camassa-Holm equation. In the end, we refer to [TB] for a pedagogical pre- 
sentation of strategies for the proof of stabihzation of PDE control systems. 

In this paper, we study the stabihzation of the ground state for cr in a neigh- 
borhood of 0. Adapting the techniques proposed in |24j . we ensure the approximate 
stabilization of the system around 4>i^a- Note that, the whole arguments hold if we 
replace the target state by any eigenstate (f>k,(7 of the system. 

1.3. Heuristic of the proof. While trying to stabilize the ground state (/)i,cr, a 
first approach would be to consider the simple Lyapunov function 

V(*) = 1-|(VI/ I <^i,,)|2. 

Just as in the finite the dimensional case [7|, the feedback law 

= 5((x* I M (<^i,. I *)) 

where 5 denotes the imaginary part of a complex, ensures the decrease of the Lya- 
punov function. However, trying to adapt the convergence analysis, based on the 
use of the LaSalle invariance principle, the pre-compactness of the trajectories in 
constitutes a major obstacle. Note that, in order to be able to apply the LaSalle 
principle for an infinite dimensional system, one certainly needs to prove such a pre- 
compactness result. In the particular case of the infinite potential well, it even seems 
that, one can not hope such a result. Indeed, phenomenons such as the L^-mass lost 
in the high energy levels do not allow this property to hold true. 

Similarly to ^24j . the approach of this paper is to avoid the population to go 
through the very high energy levels, while trying to stabilize the system around 

As in Theorem [L2l let us consider F > 0, s > 0, e > 0,7 > 0, cr G M. First, we 
consider the case , a ^0. Let G H^,^-^{I,C) with 

ll^-olU^.^^ SjFand |(^'o,</'i,.)| >7- 

We claim that there exists N — N{T, s, e, 7) G N*, large enough, so that 

00 2 

E |(*o,0,,,)|2<|7_. (16) 

k=N+l 

Then, we consider the Lyapunov function 

JV 

V(*) = 1 - I (VP I P - (1 - 6) ^ I (* I M p. (1.7) 

Note that, this Lyapunov function depends on the constants F, s, e, 7 through the 
choice of the cut-off dimension, N. Just like [24j, it encodes two tasks: 1- it prevents 
the L^-mass lost through the high-energy eigenstates; 2- it privileges the increase of 
the population in the first eigenstate. 

When ^! solves (S) with some control u ^ a + v, we have 



where 



ai := 1 and Ofc := 1 - e for fc = 2, • • • , TV. (1.8) 
Thus, the feedback law 

N 

^3 (^afe(x'i' I M{(lyk,a I (1.9) 

where > is a positive constant, trivially ensures the decrease of the Lyapunov 
function ()1.7p . We claim that, the solution of p.H) - (|1.2p - (|1.3p with initial condition 
^'o and control u — a + v{'i>) satisfies 

limsupdisti2(*(<),Ci,^)2 ^ e. (1.10) 

t — ^+oo 

Note that, the claimed result here is much stronger than the one provided in |24j for 
the finite potential well problem. In fact, here, we claim the almost global approximate 
stabilization of the system round the eigenstate (jji^^- 

The limit ()1.10p will be proved by studying the L^{I, C)-weak limits of ^'(i) when 
t —^ +00. Namely, let (in)nGN be an increasing sequence of positive real numbers 
such that tn — + +00. Since |15'(tn)||i2(7 c) = 1, there exists ^foo & L'^{I,'C) such that, 
up to a subsequence, ^'(t„) ^co weakly in i^(/,C). Using the controllability of 
the linearized system around (which is equivalent to {4>i,c7, x(f)k,a-) 7^ for every 
k G N*), we will be able to prove that "^oo = /30i,cr, where f3 G C and > 1 — e. 
This willimply (HIlOl). 

Therefore, by weakening the stabilization property (i.e. ask approximate stabi- 
lization instead of stabilization) we avoid the compactness problem evoked at the 
begining of this section. 

Note that, the controllability of the linearized system around the trajectory 
plays a crucial role here. This is why the developed techniques for a ^ can not be 
applied, directly, to the case of cr = 0. 

Now, let us study the case a = 0. As emphasized above, the previous strategy does 
not work for the approximate stabilization of because the linearized system around 
(pi is not controllable. The idea is thus to use the above feedback design (|1.9p with a 
dynamic a = a{t) that converges to zero as i — s- +00. Formally, the convergence of ^ 
toward Ci^a-{t) rnust happen at a faster rate than that of a toward zero (see Figure [TTT|) . 

In this aim, we consider the Lyapunov function 

N 

V{^) = 1 - (1 - e) ^ I I ^,,.(^)) p - 6| (* I <^i^,(«,)) (1.11) 

k=l 

where the function \I' 1— > ct(^) is implicitly defined as below 

a{^) = eiV{^)), (1.12) 

for a slowly varying real function 6. We claim that such a function (t(^) exists. When 
^ solves (S), we have 

, N s 

f = -2«(*)Q( ^ a,, (x^ I I *) j 
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Fig. 1.1. 



where 9? denotes the real part of a complex number, {ak)i^k^N is defined by (|1.8p and 
the notation '^^''^"J''^'' means the derivative of the map a i— > (pk.a taken at the point 
a = o'(\l/). By definition of cr(^'), we have 

Thus, the feedback law u(*) := + where 

N 



with <j > 0, ensures 



fe=i 



where 



- = 1 + 20'(V(vE'))5R(5]afc(^',0,,.(*))(^^,vI;)) 

is a positive constant, when j|6''||L°c is small enough. Thus t ^ V(^'(i)) is not increas- 
ing. 

We claim that, the solution of (|l.ip - (|1.2p - (|1.3p with initial condition and 
control u = (t(5') + satisfies 

limsupdistL2(1'(t),Ci)2 <^ e. (1.13) 

t — >+oo 

Again, this will be proved by studying the C)-weak limits of 'i'(i) when t +oo. 

1.4. Structure of the article. The rest of the paper is organized as follows. 

The Section [2] is dedicated to the proof of the Theorem 11.21 when cr ^ 0. We 
derive this theorem as a consequence of a stronger result stated in Theorem 1 2. II 

This theorem and a straightforward corollary (Corollarv l2.2p . leading to the The- 
orem 11.21 in the case tr 0, will be stated in Subsection 12.11 The Subsection 12.21 is 
dedicated to some preliminary study needed for the proof of the Theorem 12 . 1 1 and the 
Corollarv l2.2l The proofs will be detailed in Subsection [ 



The Section [3] is devoted to the proof of the Theorem 11.21 , in the case a = 0. 
Again, this theorem wiU be derived as a consequence of a stronger result stated in 
Theorem O 

In Subsection 13. H we state a Proposition (Proposition [ST]) ensuring the existence 
of the impUcit function cr — o'(5'). Then, we state the Theorem 13.21 and a straightfor- 
ward corollary (Corollarv l3.3p . leading to the Theorem 11.21 in the case cr = 0. A pre- 
liminary study, in preparation of the proof of the Theorem l3.2l and the Corollary [331 
will be performed in Subsection 13.21 The proofs will be detailed in Subsection 13.31 

Finally, in Section |31 we provide some numerical simulations to check out the 
performance of the control design on a rather hard test case. 

2. Stabilization of Ci_„ with a 0. 

2.1. Main result. The main result of Section [2] is the following theorem. 
Theorem 2.1. Let N e W. There exists cr' = cr'(iV) > such that, for every 
a e {-a^,a^) - {0}, 7 G (0, 1), e > 0, and *o e § verifying 

OC 2 

E |(M/o,0fe,,|2 < and |(vl/o,0i.,)| ^7, (2.1) 

k=N+l 

the Cauchy problem with u(t) = cr + v^^NA'^it)), 

Va,NA^) ■■= - e) + 0i,.)(*,0i..>j (2.2) 

has a unique weak solution . Moreover, this solution satisfies 

liminf |(4'(t),(/.i,,)P ^ 1-e. (2.3) 

t — »-t-oo 

The Theorem 12.11 provides an almost global approximate stabilization. Indeed, any 
initial condition G § such that (5'o,0i,cr) ^ can be approximately moved to 
Ci,cr- Notice that the regularity assumption '^0 G H^^{I,C) stated in Theorem 1 1.21 is 
not necessary for this purpose. Indeed, the feedback law depends on the initial state 
through the choice of the cut-off dimension N. 

The following corollary states that the quantity N appearing in the feedback law 
may be uniform when is in a given bounded subset of H^^-^{I, C). 

Corollary 2.2. Let s>0, e>0, r>0 and 7 e (0, 1). There exists a** = 
cr** (r, s, e, 7) > and N = N{T, s, e, 7) G N* such that, for every a G (— cr** , cr**) — 
{0}, and *o G H^^^(L,C) D S verifying 

W^oWh^^^ and |(*o,0i,.)| ^7, (2-4) 

the Cauchy problem il.l\) - l[l7^) - il.3\) with u — a + Uo-,A',is(^), has a unique weak 
solution "if. Moreover, this solution satisfies 112. 3\) . 

Remark 1 . The Theorem in the case a ^ is a direct consequence of the 
previous corollary. The feedback law mentionned in the Theorem \1.2\ is explicitly given 
in the Corollarv \2.2i 

Notice that, in the particular case a ^ 0, the Corollarv \2.2i is slightly more general 
than the Theorem ] 1.2[ In fact, the assumption "^0 S H"^ C\ H}^{I ,'C)"is not needed as 
we deal with weak solutions instead of strong ones. Trivially, this solution will be a 
strong solution for G n Hq{I, C). 

In the case a = 0, this will no longer be the case : we will need solutions in 
C^(M.,L'^) (for which the assumption ^0 G H Hq{I,C) is needed, see the Proposi- 
tionUn]). 

8 



2.2. Preliminaries. This section is devoted to the preHminary results, that wiU 
be apphed in the proof of the Theorem 12.11 



2.2.1. Eigenvalues and eigenvectors of Aa-. 

Proposition 2.3. For every k e W, the eigenvalue a ^ Xk,a G R and the 
eigenstate a (f)k,rT G [H^ H i?o)(/,C) are analytic functions o/ct G M around cr = 
and the expansion Xk,a = A/c + cr'^A^,^'' + o{a^) holds with 
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24A:2 87r2fc4- 



(2.5) 



There exists a* > 0, C* > such that, for every ao,(Ti £ (— cr*,(T*) — {0}, for every 
k G N% 



{X(, 


h,aai4'k,cro) 


^0, 


(2.6) 




Afc,o-o — 




C*a^ 
^ k ' 


(2.7) 




da 




C* 


(2.8) 




da 




c*, 


(2.9) 






w 


^ C*\ao-ai\ 
k 


(2.10) 



In the previous proposition, the notation means the derivative of the map 



a I—* ^feg. taken at the point a = ao- In the same way, we will use the notation — ^ 
for the derivative of the map a i >■ Xf^jj at a — ctq. 



Proof of Provosition \2.3\ : We consider the family of self-adjoint operators Aa- = 
A — ax in the space {H"^ n Hq){I,C). In this Banach space, the operator x (as a 
multiplication operator) is relatively bounded with respect to A with relative bound 
(in the sense of [21], page 190). Therefore Aa- is a self-adjoint holomorphic family 
of type (A) (see [^T], page 375). Thus the eigenvalues and the eigenstates of A^ are 
holomorphic functions of a. 

Thanks to the Rayleigh-Schrodinger perturbation theory, we compute the first 
terms of the expansions 



Afc,, 



Afc+aAi^^-t-a^Af + 



Ok.cr 



Considering the first and second order terms of the equalities Aa(j)k,tj 
\\4'k,a\\l2 = 1, we get 



= A 



k.crfPk.cr, 



X<pk 



Al'Vfe, 



b^k\^k) 



(2.11) 



dx' 



\ ^(2) 1,(1) , \(2) , 

Afe0fe +Afe 0fc + Xk Vk, 



Taking the Hermitian product of the first equality of (|2.1ip with 0^ and applying the 



parity properties of 0^, we get A 



(1) 



0. Considering the Hermitian product of the 
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first equality of (|2.1ip witli 0j, we get 

E (2.13) 



,(1) ^ V- {x(l)j,(l)k) 
^ ^ A, - A,, 



wfiere tfie sum is taken over j G N* such that the parity of j is different from the 
parity of fc. Taking the Hermitian product of the first equality of ()2.12|) with 0^ we get 
A[,^^ = — (xc/)^^-*, Using (|2.13p and the exphcit expression of {x4>k,4>j) computed 
thanks to (jl.Sp . we get 

In order to simplify the above sum, we decompose the fraction 
in the form 



2^q^\(X-qf {X + qY-J 2^ q-^ \(X - qY ' (X + q) 



5/1 1 
+ t:^ TT-r ^ + 



27g5 (X + qY) 2»q'^\{X-qY {X + qY 



2>^q^\{X-q) (X + g), 
Inserting this relation in the sum (|2.14p and simplifying, we find 

^^k^^ ^ {^^k~ Wr^k + -k^k + W^k- I^k) ^ (^-l^) 

where 

We apply, now, the following well-known relations for the Riemann ^-function: 
These relations imply 

oo 2 1 4 



thus 



5« = <! ^ for a = 1 3 5 

*^ '0 when k is even ' ' ' 
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g2 _ \ \ — w when k is odd, 5'4 — J Is ~ when fc is odd, 
^ [ X when k is even, ^ [is when fc is even. 



Inserting this in (|2.15p . we get 

The relation (|2.6p is proved in [6l Proposition 1]. The bound (|2.7p is given in [2TJ 
Chapter 17 Example 2.14, Chapter 2 Problem 3.7]. The inequality (|2.8p is proved in 
[6| Proposition 42]. The bound (|2.9p is a consequence of (|2.8|) . Indeed, considering 
the Hermitian product in L'^{I,C) of with the equation 

and using (|2.8p together with the orthogonaHty between (j)k,ao and ^^^'J^ (which is a 



consequence of ||0/c^cr||^2 = 1), we get 



da 

2 ^. /_9,2 \ /^*\2 



da 



which gives (|2.9p . Finally, (|2.10p is a consequence of (|2.8p . □ 

Proposition 2.4. Lei TV g N*. T/iere existe cr" = cf'^{N) > smc/i that, for 
every a € (-erf, tr") - {0}, j2, ^2 G N*, and ji, fci e {1, • • • , TV}, verifying ji ^ j2 and 
ki ^ k2, 

•^fei,(T ~ ^k2,cr — ^ji,a ~ ^j2,cr (2-16) 

implies (ji, ^2) = (^1,^2)- 

Proof of Proposition Let C* be as in Proposition 12.31 and c e (— o'g,crQ) 

where 

(2.17) 



" 4VC* 

First, we prove the equality (|2.16p to be impossible when j2 7^ k2 and 

-i_ 1 

max{j2,fc2}> —■ (2.18) 

We argue by contradiction. Let us assume the existence of j2,fe G N*, € 
{1, • • • , iV}, with ji 7^ j2, ki^k2, 32 ^ such that (^1^ and ((TTB)) hold. Without 
loss of generality, we may assume that max{j2, ^2} = 32 > ■ Using (|2.7p . we get 

Xn.a~Xk2.. > ^(j1-fc|)-2CV2 

^%(J|-(J2-1)2)-2CV2 
> ^(2j2-l)-2CV2, 



Using the equality of the left hand sides of these inequalities, together with (|2.17p . we 
get 

8C*crf iV2 + 1 
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which is a contradiction. 

Therefore, it is sufficient to prove the Proposition 12.41 for j2, k2 E {!,■■■ , [{N'^ + 
l)/2]}. Moreover, it is sufficient to prove that, for every ji,fci S {I,-- - ,N} and 
J2, e {f , • ■ ■ , [{N^ + l)/2]}, with ji ^ j2, fci 7^ ^2, (ji, j2) 7^ (fci, fc2), there exists 
4,fcij2,fc2 S (O'f^S) such that, for every a e (-^.fc^ J2,fe2)' fEUHl) does 
not hold. Indeed, then, the foUowing choice of a'^{N) concludes the proof of the 
Proposition 12. 4[ 

J{N) := min{a«^ ji, fci S {1, • • • , N},j2, fc^ € {1, • • • , (iV^ + l)/2}, 

(ji,i2) 7^ {ki,k2),ji ^j2,ki 7^ fca}. 

Let ji, fci e {1, • • • , N}, j2, fc2 e {1, • • ■ , (A^^ + l)/2} be such that ji ^ ja, fci 7^ ^2, 
{jiih) 7^ (^1,^2)- We argue by contradiction. Let us assume that, for every a\ > 0, 
there exists a e (— ctJ, fi}) such that (|2.f 6p holds. Using the analyticity of both sides 
in (j2.f 6p with respect to cr, at cr = 0, this assumption implies that 

\ (2) _ X (2) ^ X (2) _ , (2) 

fcl k2 jl j2 ■ 

Using (|2.5p together with a rationality argument, we get 

■p" TS" — 77 ■ 

^1 ^2 Jl J2 ^1 ^2 Jl J2 

Since ki 7^ ^2 and ji 7^ j2 , we deduce from the previous equalities that 

'^l ^2 "'2 '^1 ^2 J2 

Therefore fci = ji and ^2 = j2 , which is a contradiction. □ 
2.2.2. Solutions of the Cauchy Problem. 

Proposition 2.5. Let a- e M, iV e N*. e > 0. For every 4*0 G §, there 
exists a unique weak solution ^ of il.l\) - n^) - il.3\) with u{t) = cr + ticr,Ar,e(^'(i)), i.e. 
* € C°(R, S)nCi(R, if(o^(/, C)), f/ie egwaZitt/ (Ilip holds in H^^^{I, C) /or every t e R 
and the equality holds in §. 

Proof of Provosition\2.5\ : Let ct G M, G N, e > 0, *o G § and T > be such 

that 

TNe^"^ < 1. (2.19) 

In order to build solutions on [0, T], we apply the Banach fixed point theorem to the 
map 

9: CO([0,T],§) -> CO([0,r],§) 

where * is the solution of (frT1) - pT2)) - (fr3)) with u{t) = a + v^^NA^i^))- 

The map O is well defined and maps C"([0, r],§) into itself. Indeed, when ^ G 
C°([0, T], §), It : 1 1— > a+Vcr^N,e{^W) is continuous and thus the Proposition ! 1.11 ensures 
the existence of a unique weak solution ^E*. Notice that the map Q takes values in 

cO([o,r],s))nci([o,T],i7-„2). 
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Let us prove that 6 is a contraction of C"([0, T], §). Let e C°([0,T],S), 
Vj ■■= Wo-,Ar,e(^j), *j := for J = 1, 2 and A -.^^i - *2- We have 

A{t) = I / e-^-^-f'-^^iwixACs) + [vi - V2)x'^2{s)]ds. 
Jo 

Thanks to (|2.2p . we have ||l~(o,t) ^ ^ for j = 1,2 and \\vi — W2|1l~(o,t) ^ 
2iV|iei-6llco([o,T],L^)- Thus 

||A(t)||i2 <; / 7V||A(s)|U2+7V||a-6llco([o,T].L2)ds. (2.20) 
Jo 

Therefore, the Gronwall Lemma imphes 

II A(t)||co([0,T],L2) ^ 116 - 6llcO([0,T],L2)iVre^'^, 

and so (|2.19p ensures that O is a contraction of the Banach space C^dO, T], §). There- 
fore, there exists a fixed point ^ e C^^O, T], §) such that 6(*) = Since 6 takes 
values in C^^O, r],S) n C^([0, T], H^"Jj^), necessarily 4' belongs to this space, thus, it 
is a weak solution of (IlII])-((ll2l)-(ir3l) on [0, T]. 

Finally, we have introduced a time T > and, for every Vto € §, we have built a 
weak solution ^' € C°([0,T],S) of (IIIl])-(IIl2])-(IIl3l) on [0,T]. Thus, for a given initial 
condition *o & §, we can apply this result on [0, T], [T, 2T], [2T, 3T] etc. This proves 
the existence and uniqueness of a global solution for the closed-loop system. □ 

Proposition 2.6. Let a > 0, N e N, e > 0, (^'o)„gN be a sequence o/ S and 

e with ll^-C^IU^ 1 be such that 

lim = strongly in H-\l, C). 

n^+oo 

Let (resp. be the weak solution of l[n\) - in^) - l[Tl\) with u{t) = cr-hwa,jv,e(*") 

(resp. with u{t) = a + Va,N,ei'^°°it)) )■ Then, for every r > 0, 

lim *"(r) = ^°°{t) strongly m H-^{I,C). 



Proof of Proposition \ 2.6\ : Let us recall that the space H ^(/, C) (dual space 
of H^{I,C) for the L^(/, C)-Hermitian product) coincides with H'^^^{I,C) and that 

V2||.||h-i = ll-ll/f-1 (because ||.||/fi = \/2|| -Hhi^^ )• We introduce C > such that, 

llx^ll^-i C||^||^-i , e H-\I,C). (2.21) 
Such a constant does exist. Indeed, for every ^ G Hq{L,C), .x^ G Hq{I,C) and 

\H\\Hi, = (^jH' + edx^ ' < ||C'||l^(i + Cp) 

where Cp is the Poincare constant on /. Thus, for E H^^{I,C), we have 

l|a;¥'||H-i(/,c) = supUxip,0;^ ^ H^iI,Q, llCllffi = i} 

sup{||^||H-i||a:$||//i;e e H^il^C), Ullni = i} 

^ (i + cp)iieiiH-i. 
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In order to simplify the notations, in this proof, we write u(\E') instead of W(T,Ar,e(^')- 
We have 



+i /o e-^^(*-")[w(*"(s)) -'y(1'°°(s))]a;^'"(s)ds 
+1/0 e-*^(*-^)i;(«'°°(s))a;[«'"(s) - ^°°{s)]ds. 



Using (1221), ||*"(s)||i2 = 1, ||*°°(s)||i2 1 and the fact that <j)k,a,x(j)k^a e H^ilX) 
for fc = 1, • • • ,N, we get 

27VCa(7V)||(*"-*°°)(s)||ff-i, (2.22) 

where Ca{N) :— sup{||0fe,(T||jji(/^c); fc € {1, • ■ ■ ,^}}- The semigroup e^'"^* preserves 
the iJ^i-norm thus, using |w(^'°°(s))| N and (P?^ . we get 

+C /„'(|a| + 2iVa(iV) + Ar)||*"(5) - *°°(s)||H-ids. 
We conclude thanks to the Gronwall Lemma. □ 

2.3. Proof of Theorem [Hi] and Corollary [2721 



Proo/ 0/ Theorem \2.1\ : Let iV e N* . Let ct* > be as in Proposition [2]3] and 
crt) crt)(7v) be as in Proposition O Let a** := min{cr*, (jf }. 

Let CT e (-CT**,CT**)-{0}, 7 e (0,1), e > 0, G S with and * be the weak 
solution of (iri])-(ir2l)-(ir3]) with u{t) = a + Va.NA^{i)) given by Proposition [231 For 
V? G L2(J,C), we define 



AT 



V.,nAv) - 1 - l(<P: 01,-)!' - (1 - e) E M?- (2.23) 



fe=2 



Since * G C^{M.,H-^^ {!,€.)) and 0^,^ G H^^^{I,C), t <^ VN^aA'^it)) is C^. Using 
(|1.1[) . integrations by parts and ai := 1, Ofe := 1 — e when fc ^ 2, we get 



V..jv..(«') --23fi(^f^iafe(-zyl,^' + w,,jv,,(^')x^',0fc,,)(*,</.fe,,)), 
= -2i;,,jv,e(^(i))2. 



Thus, i 1-^ Va-.N.ei'i'it)) is a non increasing function. There exists a G [0, VCT^Ar.c(^'o)] 
such that Vcr.jv,e(^(0) ^ when t ^ +00. Since G S and (|2.ip holds we have 

V..^,.(vI/o) = 1 - (1 - 6)Ef.l - e|(*,01,.>P 

= 1 - (1 - 6) (1 - Er=A'+i <^fc,.)p) - '/'I,.)!' 

<i-(i-e) (i-|2i)-e72 

< e, 

thus a G [0, e). 

Let {tn)neN be an increasing sequence of positive real numbers such that tn — > +00 
when n +00. Since |j\['(i„)||i2 = 1 for every n G N, there exists ^"00 G L^(/, C) 
such that, up to an extraction 

*(t„) -> *oo weakly in L^{I,C) and strongly in H^^{I,C)- 
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Let ^ be the solution of 

iff = - Va^NA£.it))xi ,xe I ,te{0, +00), 

e(t,±l/2) = 0, 

Thanks to the Proposition l2.61 for every r > 0, 5'(t„+r) ^(t) strongly in H^^{I, C) 
when n ^ +00. Thus Va,N.e{'^itn + t)) ^ '^t7,N,(i{£.{''')) when n ^ +00, because 
^<y,N,ei-) is continuous for the L^-weak topology. Therefore Vcr_N.fXi {''')) = ^- Fur- 
thermore, the relation (|2.24p holds when ^E* is replaced by ^, and thus v^ N,e{^{T)) = 
and ^ solves 

i% = , X e I , t e (0,+oo), 
e(i,±l/2) = 0, 

m = *oo. 

Therefore, we have 

00 



fc=i 

The equality Wo-,a',£(C) = 0, then, gives 



N 



E E ak{^oo,<l>j,a){x^j.a,<l>k,a){^oo,^k,a)e'^^'''--^^-'>^] =0. (2.25) 



Let W(fei,fe2) := Afcj^o- - Afc^^o- for every fci, fc2 G N* and 5 := {(fci, fe); fci e {1, • • • , A^}, 
k2 G N*,A:i ^ fc2}. Thanks to the Proposition 12.41 all the frequencies ujk for K G S 
are different. Moreover, there exists a uniform gap 6 > such that, for every Lu,Li) G 
{zLllik', K G S} with CO Cj, then — a)| ^ (5. Thus, for T > large enough, there 
exists C = C{T) > such that the Ingham inequality 



2 

dt 



holds, for every {aK)K^s e /^(5,C) (see [H Theorem 1.2.9]). The equality (P?^ 
implies, in particular, 



(^'oo,0i,a)(a;0j>,(/)i,^)(^'oo,?!'i,^) = 0,Vj ^ 2. 
Thanks to jM]), we get 



(«'oo>j,a)('J'oo,0i,.) = 0,Vj ^ 2. (2.26) 

Let us prove that 

(«'oo,<^l,a) 7^0. (2.27) 

Since ||5''^ ||i2 ^ 1, we have 

Va,NA'^oo) ^ 1 - |(*°°, - (1 - e) Er=2 -/-fc,.)!^ 

= 1 - - (1 - e)[||*°°||i. - |(vI'-,(/.i,.)|2] 

^e-e|(vl/-,0i,,)|2. 
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Moreover, Vcr,N,e{'^oo) = a < e, thus 

which gives (I2.27p . Therefore (|2.26p justifies the existence of /3 e C with \/3\ ^ 1 such 
that ^-oo = /3^i,<T- Then, e> a = VN,aA'^°°) = ^ " l/^P' thus > 1 - e. Finally, 
we have 

lim |(*(t„),0i.,)|2 = |(*oo,0i..>|' = > 1-e. 

n — ^+00 

This holds for every sequence (i„)„gN thus (|2.3p is proved. □ 

Proof of Corollary \2.2\ : Let C* , cr* > be as in Proposition 12.31 There exists 
N = N{T, s, e, 7) e N* large enough so that 

(2.28) 



Let a** = <7**{N) be as in Theorem [211 (notice that a** < <t*) and cr G (-cr**, ct**)- 
{0}. Let ^0 G H^^,^{I,C) n § verifying (|^ . In order to get the conclusion of 
Corollarv l2.2l we prove that (|2.ip holds, and we apply the Theorem 12. II Using (|2.7p , 
we get 

00 



N + l, 



' A 



iv+i--jv+r 



Thus ((2:28)) implies (EH). □ 

3. Stabilization of Ci. In all this section, the constants C* , a* are as in Propo- 
sition 12.31 



3.1. Main result. First, let us state the existence of an implicit function cr(\I>) 
that will be used in the feedback law. When X is a normed space, a G X and r > 0, 
we use the notation Bx{a,r) :— {y G X; \\y — a\\x < r}. 

Proposition 3.1. Let N eW , e> 0, and 9 e C°°(R+, [0,cr*]) be such that 

9{0)=0, eis)>0 Vs>0, WO'Wl^^^^. (3.1) 
There exists a unique a G C°°{B]^2(0. 2), [0, ||0||l°°]) such that 

<j{^^)=0iV,(^).N,em^ V^G 5^2(0,2), 

where Va-,N.<: is defined by I12.23\) . 

The proof of this proposition is done in [7j. For sake of completeness, we repeat 
it in the Appendix. The main result of this section is the following. 

Theorem 3.2. Let G N*, 7 G (0,1), e > 0, 6I G C°°(M+, [0, cr*]) verifying 
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and 

Leta e C°^(Si2(0,2),[0, ||6'||loc]) be as in Proposition[EJi For every ^oeSn{H^r\ 
H^)iI,C) with 

Er=Ar+il(*o,0fc)P<32(l?W ""'^ |(*o,0i)|>7, (3.4) 

the Cauchy problem lll.l]) - [T^} - lll.3\) . with u{t) = a{^{t)) + Wo-(^(t)),Ar,e(5'(t)) has a 
unique strong solution ip. Moreover this solution satisfies 

liminf|(*(t),0i)|2 ;^ 1-e. (3.5) 

t — ^ + 00 



The following Corollary states that the quantity N appearing in the feedback law 
may be uniform in a fixed bounded subset of iJ* for s > 0. 

Corollary 3.3. Let s>0, e>0, r>0 and 7 e (0,1). There exists N = 
iV(r,s,e,7) e N* such that, for every ^0 e Sn (iJ^ n Fo^)(/, C) with *o e Hf^.^{I,C), 

||*o||H(=„, anrf|(*o,0i)| ^7, (3.6) 

the Cauchy problem U.1]) - U.S!\) - U.3\) . with u{t) = cr(*(i)) + Wo-(^(t)),A'.e(*(0) ^".s a 
unique strong solution . Moreover this solution satisfies \3.5\) . 

Remark 2. The Theorem \l.S\ with a — is a direct consequence of the Corol- 
lary \3.3l The feedback law, evoked in Theorem is implicitly given by the Corol- 
lary \EM 

3.2. Preliminaries. 



Lemma 3.4. Let N e N* , e > and 9 satisfying fO]) . There exist C{N) > 
and C{N) > such that, for all ^1, ^2 e 5^2(0, 1), 

k(a) - 't(6)| < 3N\\e'\\L^Ui - (3.7) 

k(a)-'T(e2)Kc(iv)ii^'iiL~iia-6iiH-i, (3.8) 

k<.(6).JV.e(6) - «.(6),JV,.(6)I ^(1 + 3iVC*||0'|U.^.)||a - (3.9) 
lM«i),Jv..(6)-M?2),Jv,.(6)K^(A^)liei-6lli/-i. (3.10) 



Proof of Lemma \3.Ji\ : Since N and e are fixed, in order to simplify the notations, 
we remove them from the subscripts of this proof. We have 

k(Ci)-fT(6)| ^ ll^'l|L~|V.(«,)(ei)-V.(4,)(6)|. (3.11) 

Using 



I (Ci ) I ' - 1 (6 >fc..2 ) P = (a - 6 , <i^K., ) (fi , <i^K.^ ) 
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(3.12) 



and (pmi)) . we get 

- V,(^,)(C2)| 27VCi(7V)||a - + 2iVC*|fT(ei) - 

where Ci(A^) := max{||(y9i..cr||_f/i; ^ ^ {1: ■ • • j^liO" G [O^f*]}- Using the previous 
inequalities and p.ip . we get 

i^|ai-a2K27VCi(7V)C||0'|U||ei- 611^-1, 

which imphes dSJl) and dSH) with C{N) = SNCCiiN). 

Let us write vj instead of v^(^^.-f{^j). Using, for the term 

(2^6,0j,<T(ei))(fl''^i,'T(6)) - (2^6,'/'i,a(C2))(f2,'/'jX6)) 

the same kind of decomposition as in (|3.12p . together with (|2.10[) . we get 

l^^i - V2\ < A^lia - 6I|l^ + NC*\a{^i) - 

- f2K 27VCCi(iV)||ei - + 2AfC*|(T(Ci) - (7(6)1, 

where C is defined by (g^Il)- Thus, using dST]) and ([S^), we get ([SH) and ((XTU)l with 

CiN) := 2N[CCi{N) + C*C{N)\\0' Woo]. □ 

Proposition 3.5. Let N G N* , e > 0, 9 verifying i3.1\) and For every 

eS the Cauchy problem ilT7 ]) -ilTg |) - jTg)) with u{t) = cr{^it)) + Va{^(t)),NA^i't)) 

has a unique weak solution i.e. ^ S C°(R,§) nCi((0,+oo), iJ^^p. //, moreover 

* e n i7i)(/,C), </ien * is a strong solution i.e. ^' G C"(M,i?2 n i/^) n 

Ci((0,+cx)),L2). 

Proof of Proposition \3.5\ : The strategy is the same as in the proof of Proposi- 
tion [131 Let T > be such that 

2 

Let ^0 G §■ In order to build solutions on [0, T], we apply the Banach fixed point 
theorem to the map 

9: CO([0,T],S) ^ CO([0,r],§) 

where * is the weak solution of (frT|) - (fL2| - (fL3l) with u{t) = cr{^{t})+v„(^^(^t)},N,e{i{t)). 

The map Q is well defined and maps C*'([0,T],S) into itself, moreover, it takes 
values in C°([0, T], §) n €^{{0, T),H^q^) (see Proposition [LI]). Let us prove that 6 is 

a contraction of CO([0,T],§). Let e C°([0,r],S), Vj := w^(^^.),Ar,e(Ci), := ©(Cj), 
for j = 1, 2 and A := *i - ^1*2. We have 

A(i) = I f e-'^(*-^) [(a(a) + «i)a;A(s) + {a{^i) - a(6) + vi - V2)x^2is)]ds. 
Jo 
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Using (1231) and (1X^1) . we get 



IIAWIU2 s; J^{\\e'\\L^.+N)\\A{s)\\L2ds 

+ /o {m\S'\\L^' + N[l + 3NC*\\9'\\l^]) 11^1 - ^2\\L-ds. 
Thus, the Gronwall Lemma imphes 

l|A||co([o,T],L^) ^Ui- 6llco([o,T],L^)[l + 3(1 + iVC* ) || ]iVTe^[^+ll''ll--l . 

The choice of T and p.3p ensure that is a contraction of C"([0, r],§). Therefore, 
there exists a fixed point ^ e (^"([O, T], S) such that 6(5') = Since 6 takes values 
in C"([0, T], S)nC^([0, T],H^qj), necessarily ^ belongs to this space, thus, it is a weak 
solution of dlllll-IOll-dOl) on [0,r]. 

If, moreover, *o e (-ff^nif(})(/,C), then the map 6 takes values in C"([0, T],H^r\ 
Hq) n C^([0, T], L^) thus 5* belongs to this space and it is a strong solution. 

Since the time T does not depend on 5*0, the solution can be continued globally in 
time. We, therefore, have the existence of global solutions to the closed-loop system. 
□ 

Proposition 3.6. Let a > 0, N e N, e > 0, 9 as in ^J^, (5'o)„eN be a sequence 
o/S and G ^j^^ ||5'g°||^2 ^ 1 6e such that 

lim = stronc//?/ m iJ"^(/, C). 

Let (resp. 6e t/ie w;ea/t solution of (CUP -(Eli -(Eli wi^/i w(t) = cr(*"(t)) + 

v<.(*"(t)),jv,e(*"(0) r^esp. withu{t)^a{-9^')+v,^^^),NA'^°"it)))- Then, for every 

T>0, 

lim *"(r) = *°°(r) strong;/?/ m F-\/,C). 

n — ^+oo 

Proof of Proposition \3.6\ : The proof exactly follows that of the Proposition 12.61 
In order to simplify the notations, we write w('I') instead of WCT('i').Ar,e(5'). We have 

(*" - *°°)(i) = e-*'^*(*JJ - ^-If ) + i /q e-*^(*-^) [cr(*") - cr(*°°)]x*"ds 
+i J* e-*^(*~") [«(*") - w(5'°°)]a;5'"ds 

Jo e'*^(*"") [ct(5'°°) + w(1'°°)]a;(5'" - ^'°°)ds. 

Using jUl), (13T0| and < ||x*||l2 < 1, we get 

ll(*"-*°°)(i)llH-i ^ W^^-^^Wh-^ 

+ lo {ciNWh^ + CiN)+Ci\\0\\L^ + TV)) II*" - ^°-'\\H-ids, 

where C is given by (|2.2ip . The Gronwall lemma concludes the proof. □ 
3.3. Proof of Theorem m and Corollary [3731 

Proof of Theorem [XH ; For ip £ Bl2 (0, 2), we define 

19 



where Vcr,N,e is defined by (|2.23p . Since N and e are fixed, in order to simplify the 
notations, we omit them in the subscripts of this proof, and we write instead of 

'^<T(*),Af,e(*)- 

Let e § n n H^){I,C) and * be the strong solution of (frT|) - pr2)) - (fr3l) 
withii(t) = CT(*(t))+w^(^(t))^jv,£(*(i)) given by ProDOsition[3?5l Since * e C^{R,L^) 
and a e C°°(Bl2(0, 2)), the map t ^ V{^j{t)) is C^. We have 



N 



k,{T 



dt 



\k=l 



da 



)(*,<^fe,<T(*)) 



where ai := 1 and aj, := 1 — e for k — 2, - ■ ■ , N. Moreover, 



d 
di 



thus 



N 



l + 20'(V(7^))3?^afc(*,- 
\fc=i 

Using jMl) and ^J^, we get 

/ N 

l + 20'(V(V))5R^afc(*, 



da 



)(*,</'fe,^(*)> 



V(*) = -2t>(^')2. (3.13) 



da 



<t(*) 



)(^,0fe,.(*)) ^ i~2\\e' w Nc* >o 



thus, < i-^ V(5'(t)) is a non increasing function. There exists a E [0, V(^'o)] such that 

lirr 

Using dmni), (ESI) and ([SH), we get 



lim V(«'(t)) = a. 



|(*O,01.<T(*„)>| |(*O,01>| - |(*O,01 -'/'l,<T(*o)>l 

2(c'|ie|Uoo)" 

AT 







16(l-e/2) 






16(l-e/2) 






(l-£) 



where e := e/2. Thus, as in the proof of Theorem 12.11 V(^'o) < e, so a £ (0, e). 

Let (tn)nGN bc an increasing sequence of positive real numbers such that i„ +oo 
when n +oo. Since ||5'(i„)||i2 = 1 for every n G N, there exists vj/^ g L^{I,C) 
such that, up to an extraction 

*(t„) ^ *oo weakly in L^(/,C) and strongly in iJ"i(/,C). 

Let ^ be the weak solution of 

e(t,±l/2) = 0, 

m = *oo. 
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Thanks to the Proposition l3.61 for every t > 0, 5'(t„+T) ^(t) strongly in H ^{I, C) 
when n +00, thus o'(5'(t„ + r)) when n +00 (see Lemma I3.4p . 

Therefore, V(4'(in + r)) ^ V(C(t)) when n +00, so V(0 = a. Thus, cr(^) ee := 
6'(a) and we have, for every t £ R+, 



AT 

2 



Vim) = 1 - K^W, - (1 - E K^W' '^^-^ 



fc=2 



Since ^ G (R+ , -ff^p^ , the previous equahty imphes 



dt 



Since V(C) = a, then v{Cj = 0. 

Firsf case ; a = 0. Then V(5'(t)) — > when t +00 and a — 0. Moreover, for 
every t G (0, 00), 

v(*(t)) ^ 1 - </>i,.(*))p - (1 - e) Er=2 

^e(l-|(*,(/.i,,(*))|2). 

Thus, 

|(^'(i„)>i,^(*(t„)))l ^ 1, 

which leads to 



|(vl/(i„),0i)|->l- 

because (T(4'(t„)) ^ 0. 

Second case : a ^ 0. Then a — 6{a) > 0. Exactly as in the first analysis, done 
in the proof of Theorem 12. ![ we get 

where /3 e C and \f3\'^ > 1 - e. Thus 

lim |(^'(i„),0i)| = \{^oo,<l)i)\ > l/SI-K^'oo^i,?-^!)! > Vl-e/2-C*a 

n — *+oo 

where we used p.7p in the last inequality. Finally, thanks to <a^ \\(^\\oo and (|3.2p . 
we get (|3?5|) . □ 

Proof of Corollary \3.3\ : It can be done in a very similar way to the proof of the 
Corollary O □ 

4. Numerical simulations. In this section, we check out the performance of 
the techniques on some numerical simulations. We consider, as a test case, the sta- 
bilization of the initial state ^'o — ■^{(j)i,a + (f>3.<j) around the ground state (j)!,^- 
Therefore, the cut-off dimension is 3. Note that, such a test case is particularly a 
hard one in a near-degenerate situation. Indeed, considering the feedback law ()1.9p 
for (7 = 0, one can easily see that for parity reasons = 0. 

In a first simulation, we consider the non-degenerate case of a ^ 0. As mentioned 
above the constant a needs to be small. In fact, one should choose a, such that the 
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Fig. 4.1. The approximate stabilization o/Ci,ct, where 'i/o = -^((pi^a +<t>3,a) and therefore the 
cut-off dimension is 3; as it can be seen the closed-loop system reaches the . 05-neighborhood of ij>i,a 
in a time T = ISOtt corresponding to about 200 periods of the longest natural period corresponding 
to the ground to the first excited state. 

perturbation ax is small compared to the the operator — ^ . We choose it here to 
be cr = 2e + 01. The Figure |4?T] illustrates the simulation of the closed- loop system 
when u = a + with <; = le + 03 and e — be — 02. The simulations have been 
done applying a third order split-operator method (see e.g. |T6j), where instead of 
computing exp(— i dt (Aa — v^x)) at each time step, we compute 

exp(— i dt A(j /2) exp(i dt v^x) exp(— i dt A^/2). 

Moreover, we consider a Galerkin discretization over the first 20 modes of the sys- 
tem (it turns out, by considering higher modal approximations, that 20 modes are 
completely sufficient to get a trustable result). 

Now, let us consider the degenerate case of ct = 0. As mentioned above, such a case 
is not treatable with the explicit feedback design of (|1.9[) . However, the simulations of 
Figure l472j show that the implicit Lyapunov design provided in Subsection ll.3l removes 
the degeneracy problem and ensures the approximate stabilization of the initial state 
-^(01 + 4>3) around the ground state 4>i. 

We consider the function 9(r) = rjr with 77 = 7e-|-02. Furthermore, in the feedback 
design Ve, we consider ^ = le -I- 03 and e = 5e — 02. The numerical scheme is similar to 
the simulations of Figure 14.11 In order to calculate the implicit part of the feedback 
design (t(^'), we apply a fixed point algorithm. 

5. Appendix. This appendix is devoted to the proofs of the Proposition ! 1 . II and 
the Proposition [2t1 

5.1. Proof of Proposition HH] . Let *o e §, Ti > and u e C°{[Q,Ti],R). 
Let T e (0, Ti) be such that 



II"I|l1(0,T) < 1- 
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Fig. 4.2. The approximate stabilization of Ci, where >I'o = ~^i4'i ~^ <t>3) "■'"■^ therefore the 
cut-off dimension is 3; as it can be seen, the closed-loop system reaches the . 05-neighborhood of ij>i 
in a time T = IOOOtt corresponding to about 1300 periods of the longest natural period corresponding 
to the ground to the first excited state. 



We prove the existence of * e C''([0, T], C)) such that (fL4ll holds by applying 
the Banach fixed point theorem to the map 

9: C7"([0,T],L2) ^ CO([0,r],L2) 
where ^ is the weak solution of 

^{t,±l/2) = 0. 
i.e. 5- e C°([0,T],L2) and satisfies, for every t £ [0,T], 

4-(t) = e-'^'^Q + I f e-''^'^'-''>u{s)x£,{s)ds in L^I, C). 
Jo 

Notice that 6 takes values in C^{[0, T],H^^^{I, C)). 

For a, 6 e C0([0,r],L2(/,c)), *i :=e(a), *2 := 0(6) we have 

(*! - ^2){t) = I f e-^^(*-^'u(s)x(Ci - ^2){s)ds 
Jo 

thus 

|l(*l~*2)(t)||L^ / Hs)\dsUi-^2\\c»ilO.T].L-)- 

Jo 

The assumption (jS.ip guarantees that 8 is a contraction of C°([0, T], L^), thus, O 
has a fixed point * e C"([0, T], L^). Since 9 takes values in C^{[0,T], H^^^^), then ^ 
belongs to this space. Moreover, this function satisfies (|1.4p . 
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Finally, we have built weak solutions on [0, T] for every ^fp, and the time T does 
not depend on 5*0, thus, this gives solutions on [0, Ti]. 

Let us prove that this solution is continuous with respect to the the initial con- 
dition VPo, for the L^(/, C)-topology. Let ^'o,^'o G § and <f> the associated weak 
solutions. We have 

thus Gronwall Lemma gives 

IK* - <^){t)\\L2 ^ 11*0 - $o||L2e"""^^<°.^i). 

This gives the continuity of the weak solutions with respect to the initial conditions. 

Now, let us assume that *o G H Hq{I, C). Take C to be a positive constant 
such that for every Lp e H'^r\Hl{I, C), \\x(p\\u2f^ui ^ C'|l</'|l//2n^i . We consider, then, 
T > such that C||u||ii(o,T) < 1- By applying the fixed point theorem on 

92 : C\[0,T],H^r^Hl) ^ C°([0, T], i?^ n iJ^) 

defined by the same expression as O, and using the uniqueness of the fixed point of 
0, we get that the weak solution is a strong solution. The continuity with respect 
to the initial condition of the strong solution can also be proved applying the same 
arguments as in above. 

Finally, let us justify that the weak solutions take values in S. For G H 
Hi, the solution belongs to C^{[0,T],L'^) n C°([0,r],H2 n Hi) thus, the following 
computations are justified 

Thus e § for every t e [0, T]. 

For *o G we get the same conclusion thanks to a density argument and the 
continuity for the C''([0, T], L^)-topology of the weak solutions with respect to the 
initial condition. □ 

5.2. Proof of Proposition [37l1 . Let * e (0, 2). We prove the existence of 
(t(*) by applying the Banach fixed point Theorem to the map 

n: [0,||^|U»] ^ [OM\l^] 

For cri,o-2 G [0, we have 

|n(ai) - n(a2)| ||0'||l»|V.„jv,.(*) - V.,.jv..(«')|. 
Using the following inequality 

|K*,0J,<^i>|' - \{'^,4>J.<ro_}\'\ S;|(*,9ij,<.i -0,,a2>(*,0J,<.i>| +\{'^,4>J.<r2)T^Aj77~K^\ 

^8\\<t>3.Ti - 4>j,:72\\l2, 
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together with (|2.10p . we get 



Thus, the assumption (|3.ip ensures that 11 is a contraction of [0, ||0||ioc]. Therefore, 
n has a unique fixed point cr(^'). 

Now, let us prove that a is . The map 

F: [O,||0||loo] X 5^.(0,2) ^ R 

is regular with respect to a and F((7(^'), ^I*) = 0, for every 5* e Bjjii^, 2), and 

1 



— *) = f - 20'(V,(*),Ar,,(*))— \^a,N,m 

Indeed, for ctq G [0, |16'||l~] and * G 5^2(0,2), we have 



T(V) 2 



(5.2) 



N 



fe=i 



da 



»fe,CT0, 



where oi := 1 and := 1 — e for fc = 2, .., N . Thus, using (|2.8|) . we get 

d r. 



(9cr 



s$ 87VC* 



We get the inequality in (|5.2p thanks to the previous inequality and p.ip . 

For every e i3i2(0, 2), the implicit function theorem provides the existence of a 
local C°° parameterization cr(^) for the solutions of F((j(^^),i^ = 0, in a neighborhood 
of 4". The uniqueness of the fixed point (t(^) justifies that a and a coincide, thus a is 
C°°. □ 

Acknowledgments : The authors thank J-M. Coron, R. van Handel, O. Kavian 
and P. Rouchon for helpful discussions. 



REFERENCES 

[1] F. Albertini and D. D'Alessandro. Notions of controllability for bilinear multilevel quantum 
systems. IEEE Transactions on Automatic Control, 48(8):1399-1403, 2003. 

[2] C. Altafini. Controllability of quantum mechanical systems by root space decomposition of 
su(n). Journal of Mathematical Physics, 43(5):2051-2062, 2002. 

[3] A. Elgart and J. E. Avron. Adiabatic theorem without a gap condition. Commun. Math. Phys., 
203:445-463, 1999. 

[4] L. Baudouin, O. Kavian and J. -P. Fuel. Regularity for a Schrodinger equation with singular 
potentials and application to bilinear optimal control. Journal of Differential Equations , 
216 Issue 1: 188-222, 2005. 

[5] L. Baudouin and J. Salomon. Constructive solutions of a bilinear control problem. C. R. Math. 
Acad. Sci. Paris , 342 (2): 119-124, 2006. 

[6] K. Beauchard. Local controllability of a 1-D Schrodinger equation. Journal de Mathematique 
Pures et Appliquees, 84:851-956, 2005. 

[7] K. Beauchard, J.-M. Coron, M. Mirrahimi, and P. Rouchon. Implicit Lyapunov control of finite 
dimensional Schrodinger equations. System and Control Letters, 56: 388—395, 2007. 

[8] K. Beauchard and J.M. Coron. Controllability of a quantum particle in a moving potential 
well. J. Functional Analysis, 232: 328-389, 2006. 

[9] Y. Chen, P. Gross, V. Ramakrishna, H. Rabitz, and K. Mease. Competitive tracking of molec- 
ular objectives described by quantum mechanics. J. Chem. Phys., 102:8001-8010, 1995. 



25 



[10] C. Cohen- Tannoudji, B. Diu, and F. Laloe. Mecanique Quantique, volume I&; II. Hermann, 
Paris, 1977. 

[11] J.M. Coron. Global stabilization for controllable systems without drift. Math. Control Signals 

Systems, 5:295-312, 1992. 
[12] J.M. Coron. On the null asymptotic stabilization of the 2-D incompressible Euler equation in a 

simply connected domain. SIAM J. Control and Optimization, 37 (6), 1999, p. 1874-1896. 
[13] J.M. Coron Control and nonlinearity. Mathematical Surveys and Monographs. 136, 2007. 
[14] J.-M. Coron and B. d'Andrea-Novel. Stabilization of a rotating body-beam without damping. 

IEEE Transactions on Automatic Control, 43 (5), 1998, p. 608-618. 
[15] J.-M. Coron, B. d'Andra-Novel and G. Bastin. A strict Lyapunov function for boundary control 

of hyperbolic systems of conservation laws. IEEE Transactions on Automatic Control, 52 

(1), 2007, p. 2-11. 

[16] M. D. Feit, J. A. Fleck, Jr., and A. Steiger. Solution of the Schrodinger equation by a spectral 
method J. Comput. Phys., 47, 412-433, 1982. 

[17] O. Glass. Asymptotic stabilizability by stationary feedback of the two-dimensional Euler equa- 
tion: the multiconnected case, SIAM J. Control Optim. AA (2005), no. 3, pp. 1105-1147. 

[18] O. Glass. Controllability and asymptotic stabilization of the Camassa- Holm equation, preprint 
2007. 

[19] R. Van Handel, J.K. Stockton, and H. Mabuchi. Modeling and feedback control design for 
quantum state preparation. J. Opt. B: Quant. Semiclass. Opt. (special issue on quantum 
control), 7:S179-S197, 2005. 

[20] S. Haroche. Controle de la decoherence: theorie et experiences, 2004. Notes de cours. College 
de France, http:/ / www.lkb.ens.fr / recherche / qcdcav/collcgc/colle ge.htmI| 

[21] T. Kato. Perturbation Theory for Linear Operators. Springer, 1966. 

[22] W. Krabs On moment theory and controllability of one dimensional vibrating systems and 

heating processes. Springer Verlag, 1992. 
[23] B. Li, G. Turinici, V. Ramakrishna, and H. Rabitz. Optimal dynamic discrimination of similar 

molecules through quantum learning control. J. Phys. Chem. B., 106(33):8125— 8131, 2002. 
[24] M. Mirrahimi. Lyapunov control of a particle in a finite quantum potential well. In CDC06, 

San Diego, 2006. 

[25] M. Mirrahimi and R. Van Handel. Stabilizing feedback controls for quantum systems. SIAM 

Journal on Control and Optimization, 46 (2): 445-467, 2007. 
[26] M. Mirrahimi, P. Rouchon, and G. Turinici. Lyapunov control of bilinear Schrodinger equations. 

Automatica, 41:1987-1994, 2005. 
[27] M. Mirrahimi, G. Turinici, and P. Rouchon. Reference trajectory tracking for locally designed 

coherent quantum controls. J. of Physical Chemistry A, 109:2631—2637, 2005. 
[28] A. Pazy. Semigroups of Linear Operators and Applications to Partial Differential Equations, 

volume 44 of Appl. Math. Sci. Springer, 1983. 
[29] V. Ramakrishna, M. Salapaka, M. Dahleh, and H. Rabitz. Controllability of molecular systems. 

Phys. Rev. A, 51(2):960-966, 1995. 
[30] P. Rouchon. Control of a quantum particle in a moving box. Technical Report A/325, Ecole 

des Mines de Paris, Centre Automatique et Systemes, October 2002. 
[31] S. Shi, A. Woody, and H. Rabitz. Optimal control of selective vibrational excitation in harmonic 

linear chain molecules. J. Chem. Phys., 88(11) :6870-6883, 1988. 
[32] M. Sugawara. General formulation of locally designed coherent control theory for quantum 

systems. J. Chem. Phys., 118(15):6784-6800, 2003. 
[33] H.J. Sussmann and V. Jurdjevic. Controllability of nonlinear systems. J. Differential Equations, 

12:95-116, 1972. 

[34] G. Turinici. Controllable quantities for bilinear quantum systems. In Proceedings of the 39th 
IEEE Conference on Decision and Control, pages 1364—1369, 2000. 

[35] G. Turinici and H. Rabitz. Wavefunction controllability in quantum systems. J. Phys. A, 
36:2565-2576, 2003. 



26 



